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Abstract. Minia ti & Bell ( 2011| ) proposed a mechanism for the generation of mag- 



netic seeds that is based the finite resistivity of the low temperature IGM in the high 
redshift universe. In this model, cosmic-ray protons generated by the first generation of 
galaxies, escape into the intergalactic medium carrying an electric current that induces 
return currents, j t , and associated electric fields, E = rjj t there. Because the resistivity, 
7], depends on the IGM temperature, which is highly inhomogeneous due to adiabatic 
contraction and shocks produced by structure formation, a non-vanishing curl of the 
electric field exists which sustains the growth of magnetic field. In this contribution we 
have developed an approximate numerical model for this process by implementing the 
source terms of the resistive mechanism in the cosmological code CHARM. Our numer- 
ical estimates substantiate the earlier analysis in |Miniati & Bell| ( |2011| ) which found 
magnetic seeds between 10" 18 and 10" 16 Gauss throughout cosmic space at redshift 
z ~ 6, consistent with conservative estimates of magnetic fields in voids at z ~ from 
recent gamma-ray experiments. 



1. Introduction 



Magnetic fields are an intrinsic part of astrophysical plasmas, though in many respects 
their origin and role is still subject of intense investigation. Recent studies combining 
gamma-ray measurements at different energies indicate the existence of magnetic fields 
in cosmic voids ( |Neronov & Vovk||2010||Taylor et al.||2011||Tavecchio et al.||2011| ). If 
confirmed (see e.g. Broderick et al.| ( |2011| )), the origin of such magnetic fields poses 
a non trivial challenge for astrophysical scenarios: the mechanisms connected to the 
Weibel instability ( |Schlickeiser & Shukla||2003[ |Medvedev et al.||2006| ) and the Bier- 
mann's battery at shocks ( Kulsrudet al. 1997 ) generate magnetic fields inside collapsed 
structures, such as filaments and clusters, not in the voids. The Biermann's battery acti- 
vated by ionization fronts ( [Gnedin et al.|2000| ) generates fields that are too small in the 
low density regions. Galactic outflows can in principle export high magnetic flux (e.g 
Bertone et al.|2006[ ), but despite valuable efforts, the hydrodynamic transport that deter- 
mines its filling factor in the intergalactic space remains highly uncertain ( |Rees|2006| ). 



Miniati & Bell ( 201 1|) proposed a mechanism to generate widespread intergalac- 



tic magnetic fields that is associated with escape of high energy particles (heretofore 
CR for cosmic-ray) from the first generation of galaxies, also responsible for the re- 
ionization of the intergalactic medium (IGM). Cosmic re-ionization is mostly driven 
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by ionizing photons produced by massive OB stars. When these stars come to an end, 
they explode as supernova whose blast wave accelerates copious CR ( |Krymskyl|1977 



Axford et al.[T977||Bell|1978t|Blandford & Ostriker|1978| ). Owing to their much larger 



diffusive mean free path compared to thermal particles, CR eventually escape from the 
parent galaxies into the IGM ( |Schlickeiser||2002| ). The charged CR protons carry an 
electric current in the IGM, j c , which is not balanced by the CR electrons, typically 
fewer in number by two orders of magnitude ( Schlickeiser 2002) and subject to mush 



more severe energy losses. Due to charge imbalance and induction effect a return cur- 
rent, j u carried by the cold thermal plasma is generated, which tends to cancel the CR 
current, i.e. j t ^ -j c . Due to the finite resistivity of the IGM plasma, an electric field, 
E = r/j t , is required to draw the return current. This field opposes the propagation of the 
CR but is too weak to prevent their escape. The plasma resistivity in an unmagnetized 
medium is determined by its collisional rate which only depends on temperature, and 
can be written as rj(T) = t]o(T/K)~ 3 ^ 2 (Spitzer 1956), where 770 is a constant. Owing to 
the formation of structure in the universe, driven by gravitational instability, consider- 
able density and temperature inhomogeneities existed in the IGM long before the onset 
of re-ionization. Thus the electric field is also inhomogeneous and in particular 

3 VT 

VxE = Vx7 1 j t * --T]j t x — * 0! (1) 

where we have neglected the term oc V x j t and have used Spitzer' s expression for 
the resistivity. The important point is that, since the current is unrelated to the IGM 
inhomogeneities, the curl of the electric field is non vanishing. The rotational compo- 
nent of the electric field sustains Faraday's induction and generates magnetic field. The 



analysis of Miniati & Bell] ( |2011| ) used the observed UV luminosity function of high 



redshift galaxies to estimate the production rate of CR in those galaxies and the ensu- 
ing return currents in the IGM, j t . They also used simulations of structure formation to 
estimate the temperature gradient scale-length, T/VT, entering the electric field curl. 
They found that magnetic field is robustly generated throughout intergalactic space at 
rate of 10" 17 - 10" 16 Gauss/Gyr, until the temperature of the intergalactic medium is 
raised by cosmic re-ionization when the age of the universe is roughly t =lGyr. This 
value is consistent with the conservative values implied by gamma-ray measurements. 

In this contribution, we present results from a cosmological numerical model that 
also includes the source terms associated with the resistive mechanism described above. 
The numerical model and results are described below but we anticipate that they val- 



idates the quantitative findings of Miniati & Bell (2011 ). The rest of this manuscript 



is organized as follows. The numerical model is described in Sec. |2| the results are 
reported in Sec. [3] and Sec.|4]contains summary and conclusions. 



2. Approximate Numerical Model of the Resistive Mechanism 

We start from a numerical model of cosmological structure formation. This includes the 
evolution, in an expanding background space-time, of a collision-less dark matter fluid 
and a collisional baryonic gas, coupled through self-gravity, and subject to appropriate 
initial conditions. We then extend the hydrodynamics to include the effects of magnetic 
field in the MHD approximation and modify the RHS of the governing equations to 
include the various effects due to the return current j t (see below). Because j t ^ -j c , 



for convenience in the following we express the source terms associated to the return 
current in terms of -j c . 



The full set of equations solved in a cosmological model is given, e.g., in |Miniati & 
Colella| ( [2007| ). Here we only describe the equations of hydrodynamics, which contain 
the modifications required by the resistive model. As usual, these equations employ 
comoving density, p, and pressure, p g , and peculiar velocity, u. If a(t) is the spatial 
scale factor of the universe and a its rate of change, density and pressure are modified 
as p <— a 3 p, p g <— a 3 P respectively, to scale out the effect of Hubble expansion. 
Likewise the peculiar velocity does not include the Hubble flow component (see, e.g., 
Miniati & Colella|2007[ for details). For numerical reasons we find most convenient to 
rescale the magnetic field as, B <— a^Ji, j c <— #2j c , even though it does not scale out 
completely the adiabatic expansion effects. So finally the governing equations read 



JL + -V 
ot a 



(pu) = 0, 



dpu 1 

ot a 

dpe 1 

ot a 



pu : u + p 
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u(pe + p)-— B(B-u) 
4n 



= -2-pe - -pu • - u • (- X B I + 777^.(4) 
a a \ c 



Here p = p g + ^B • B is the total sum of the gas and magnetic pressures, e = 



u • u + 



B is the total specific energy density. The thermal energy is related to the 
pressure through a y-law equation of state, e t h = Pglpij - 1), is the gravitational 
potential, -j c the return current and r\ the resistivity. The last term on the RHS of 
the momentum equation describes the interaction between the magnetic field and the 
return current, which also contributes the term before the last in the energy equation. 
The last term in the energy equation represents ohmic dissipation of the return current. 
The electric field, //jc, driving this current does not appear in the momentum equation 
as it is balanced by frictional forces, the same ones that are in fact responsible for the 
dissipation. 

The magnetic field evolution is described by Faraday's equation, with the electric 
field given by Ohm's law. The only electric fields are those induced by motions of the 
magnetized fluid as well as the return current, i.e. E = -(u/c)xB-7/j c . So the modified 
induction equation (again in comoving quantities) reads 

m i. . i. . 



at 



= -V x (u x B) + -cV x 77j c 

a a 



2 a 



The full set of cosmological equations is solved with the code CHARM (Mi niati & Colella 
2007| ). The MHD solver is based on the constrained- transport scheme and is described 
in details in [Miniati & Martin| ( |2011| ). The source terms associated with the resistive 
process have been implemented using a second order accurate time-centered scheme. 

Concerning the CR current, j c , as discussed in |Miniati & Bell| ( [201 it is typically 
dominated by the most luminous nearb y galaxy. For a bright L* ga laxy at redshift z - 6, 
with star formation rate of 6.5 M Q yr _1 Bouwens et al. (2007 2010), assuming a Salpeter 



initial-mass-function and a 30% conversion efficiency of supernova energy into CR, it 
is found 
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Figure 1 . Simulation snapshots for the distribution of baryonic gas (left) and mag- 
netic field strength (right) on a slice through the simulated volume at redshift z = 10. 
Gas density is in units of the average total (baryonic plus dark matter) matter density, 
and magnetic field is in Gauss. The physical length each panel's side is 12.7 kpc. 



where d is the physical distance from the parent galaxy and p m [ n is the minimum en- 
ergy of the escaping CR. In the following we consider a distance d ~ IMpc, which 
typically separates bright V galaxies at z ~ 6, and a volume of linear size L « so 
that the distance dependence of the CR current can be neglected within the simulated 
volume. In addition, the spatial inhomogeneities of j c are neglected, so j c is effectively 
a constant, whose value is set in accordance with Eq. [6] We neglect the effects due 
to the generated magnetic field on the CR propagation. Though they introduce some 
degree of diffusivity in the CR propagation, the net effect on the current should remain 
negligible for a field strength below 10" 15 G. 



3. Results 

Using the model described in the previous section we have run a model of cosmic struc- 
ture formation that includes the resistive generation of magnetic field. For the purpose, 
we have adopted a concordance A-CDM based with cosmological parameters given by 
the results of the WMAP experiment |Komatsu et al.| p009). We use a computational 
box with 512 3 resolution elements and as many collision-less particles to describe the 
dark matter component. For the physical size of the computational domain we consider 
two cases. In the first one the simulated volume corresponds to a cube of side L = 140 
comoving kpc. At redshift z ~ 10, corresponding to the epoch of the results presented 
below, this translates into a physical size of 12.7 kpc and a spatial resolution of 25 pc, 
sufficient to resolve temperature structures on physical scales of a few kpc, identified 
in |Miniati & Bell| ( |2011| ) as the most relevant for the resistive process. Fig. [T] shows 



a snapshot of the distribution of the baryonic gas in units of the total matter density 
(left) and magnetic field strength in Gauss (right), on a two-dimensional slice passing 



Figure 2. Histrogram of magnetic field distribution as a function of its strength at 
z = 10 (left) and z = 6 (right). 



through the simulated volume at redshift z ~ 10. As anticipated above, the matter den- 
sity distribution shows considerable structure, which translates into a corresponding 
temperature structure due either to adiabatic compression or shocks. Notice that the 
magnetic field ranges from 10" 18 to several times 10" 17 Gauss. An important feature 
is that the magnetic field with strength in this range fills is generated throughout the 
simulated volume. This is better illustrated in the histogram on the left hand side of 
Fig. [2| which shows the distribution of magnetic field as a function of magnetic field 
strength. According to the histogram, all space is magnetized and the average magnetic 
field value is of order 10" 17 Gauss. 

Given the small box size, at later time the formation of structure quickly saturates 
and the numerical model becomes inadequate to describe the intergalactic medium. 
We have thus run a second model which employs a box ten times as large as for the 
previous case, L = 1.4 comoving Mpc. We run this model down to redshift z ~ 6, when 
cosmic re-ionization raises the IGM temperature, effectively shutting down the resistive 
mechanism. Since we employ the same numerical resources as before, so the resolution 
is at least 390 kpc physical, which is still sufficient to resolve the desired temperature 
structures. The histogram on right hand side of Fig. [2] shows the distribution of the 
magnetic field strength in this larger simulated volume at redshift z ~ 6. It shows that 
the magnetic field continues to grow in strength from redshift 10 to 6 and at this later 
time it fills the intergalactic volume with a strength between 10" 18 and 10" 16 . The 
histogram peaks at a value of several times 10" 17 , which is quantitative consistent with 
the estimates in |Miniati & Bell| ( |20TT] ). 

Finally, although the governing equations l2p] include the ingredients necessary 



to drive non resonant amplification of magnetic field (BelT| |2004 2005), there is no 
appreciable effect is our simulations. This is because for the tiny resistive magnetic 
field the growth timescales of the non-resonant mechanism is too slow on the spatial 
scales resolved by our simulation. It remains unclear whether or not this mechanism 
may efficiently kick in at later times. 
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4. Summary and Conclusions 



Minia ti & Bell| ( |20lT| ) proposed a mechanism for the generation of magnetic seeds that 



is based the finite resistivity of the low temperature IGM in the high redshift universe. 
In that model, CR escaping from the first generation of galaxies, induce a return current 
and an electric field in the IGM plasma. The electric field, E = rj(T)J, depends on 
the IGM temperature, which is highly inhomogeneous due to adiabatic contraction and 
shocks produced by structure formation. It is easy to show that this electric field has 
a non vanishing curl which sustains the growth of magnetic field. The analytic results 
in |Miniati & Bell| ( [201 1[ ) were consistent with conservative values of magnetic fields in 
voids inferred cosmic-voids by gamma-ray experiments. Here we have implemented 
the source terms describing the resistive process in the cosmological code CHARM. To 
keep the problem tractable we made a number of reasonable approximations. Neverthe- 
less, the numerical model allows a more self-consistent estimate of the magnetic field 
generation by the resistive process. We have carried out numerical experiments using 
two different computational boxes to best estimate the generated magnetic fields at dif- 
ferent cosmic epochs of interest. In conclusion the numerical estimates substantiate the 
earlier analysis in |Miniati & Bell| ( |2011| ), showing that most of the made approxima- 



tions and assumptions were essentially correct. 
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